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We examine quantum decay of localized vibrations in anharmonic crystal lattice. The theory 
which describes two-phonon anharmonic relaxation can be applied both to local modes associated 
with substitutional impurity and to intrinsic local modes (ILM) in perfect lattices. It is found that 
for sufficiently high initial excitations relaxation of vibrations is non-exponential, it demonstrates 
explosion-like behavior at specific stages of evolution. The course of the relaxation is determined 
by the initial value of energy, temperature, direction of vibrations. As an example we present the 
^ ' results of calculations of the relaxation of an odd local (impurity) mode in a simple cubic lattice 

and discuss the influence of quantum fluctuations on the stability of the ILM in one-dimensional 
j /^ ■ monatomic chain. 

t> ; 

, 1 I. INTRODUCTION 

>\ 

The localization and energy transport in nonlinear systems have received much attention recently. Numerical 
techniques play an essential role in the investigation of anharmonic discrete lattices, allowing the variation of model 
parameters and the observation of different properties. One ofptlie effects - the existence of stationary localized 
vibrations in a perfect lattice (s.c. intrinsic local modes (ILM))EJ~El is of particular interest as it links local lattice 
' dynamics with physics of solitons and underlines in general the importance of strongly excited local modes. 
^\ , So far the research of strong anharmonic effects in lattice dynamics as well as molecular dynamics simulations of the 
■ kinetics of the decay of vibrational excitations^ has been carried out within the frame of classical mechanics. However, 
the account of quantum effects may be important, as it leads to new channels in the decay. In this context the problem 
of strong local vibration is also of significant interest, it gives an example when anharmonic and quantum effects can 
be desc.riJb.ed analitycally with the application of nonperturbative theory. One of the results of the corresponding 
theoryBillj is an explosion-like release of energy by the local vibration with the emission of a burst of phonons. This 
release takes place if the energy of the local mode reaches some definite critical values. Such a behaviour is new for 
£j , vibrational systems andii. is a consequence of a mutual interplay between anharmonic and quantum effects. 

The method used inoE3 is based on the assumption that a local mode is strongly excited initially and it can be 
. £h ! considered classically. Phonons are supposed to be not excited (or to be in thermal cquilibrum) . According to 
classical mechanics, the state with non-excited phonons is stable, since for the phonon coordinates qi at rest = 0) 
the anharmonic interaction is turned off. Quantum fluctuations, coming from zero-point vibrations, turn on the 
interaction and result the relaxation of the mode. .^However, the standard quantum theory, based on the description 
of anharmonic interaction as a small perturbationtl, cannot be used for the obtaining of the solution to the problem 
of relaxation, because the interaction is not weak (due to the strong excitation of the local mode). 
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II. GENERAL THEORY 



Let us consider a local vibration anharmonically interacting with phonons. The peculiarity of the problem is that 
the strong local vibration causes periodic time-dependence of the local force constants and therefore of the zero-point 
energy of the phonon system. This time-dependence results the generation of phqepns. The mechanism of this 
process has an analogy with the black hole emission mechanism proposed by Hawkingjij and with Unruh radiationlij. 
In all these cases the time dependence of the zero-point energy causes the transformation of the initial creation and 
destruction operators in time. Namely, according to Hawkingii^ 

b i =^2(a ij a j - fiijof), (1) 
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where a and b are initial (incoming) and final (outgoing) time-dependent destruction operators. It means that when 
there are no incoming particles the number operator of the i-th outgoing state is 



JVi=<Q|S+6i|0>=X; I Ail 2 , (2) 

3 

i.e. the number of particles created and emitted in a gravitational collapse can be determined calculating the co- 
efficients Pij. There exist several situations in quantum field theory where the phenomena similar to black hole 
evaporation appear. The most famous are: 



• Pair production in a static electric field revealed by Heisenbergtil. Vacuum of quantum field theory is unstable 
against the creation of charged pairs. 

• Accelerated systems which become spontaneously excited in Minkowski-space and accelerated mirrors. 

It is also known that the rate of emission of the black hole is very high near the end of its life. As a black hole emits 
radiation it loses mass. This in turn increase the rate of emission. Near the end of its life about 10 30 erg would be 
released in 0.1 s. So the emitting rate of the black hole grows in time and at some critical moment a large amount of 
energy will be released in very short time (explosion). We want to point out the analogy of this effect in the system 
of phonons in nonlinear crystal lattice. 

Cubic anharmonic interaction of the local mode with crystalline is the sum of terms ~ Xi Xy Q(t), where Xi and Xy 
are the operators of normal coordinates of crystalline phonons, Q(t) is the time-dependent classical amplitude of the 
local mode. This interaction leads to the following time-dependent Hamiltonian of the phonon system: 

H P h(t) = - ^2(xi + ijj} x\) + -Q(t) 22(eiWev)xiXii . (3) 

i ii' 

Here is the frequency of the mode i, fewey) = J2 n e in&VnWnm n = 1, 2, . . . , no, J2i Zin^m' = <W, n is the number 
of configurational coordinates contributing to anharmonic interaction. 
The Hamiltonian © can be diagonalized as followa3: 

H ph (t) = ]T nn,(t) (bf(t)bj(t) + i) , (4) 

3 

where Qj(t) are time-dependent phonon frequencies, and 

h (*) = E + Vi i ^ & t) > ( 5 ) 

are time-dependent destruction operators; expressions for Clj, /j,^ and are given inEl. 

The anharmonic interaction considered causes not only time dependence of phonon frequencies but also changes 
of phonon operators in time. The relation (fJ) is analogous to the relation of Hawking ([j]) between field operators in 
different times of gravitationally collapsing star. This is why the mechanism of a local mode relaxation is analogous 
to that of black hole emission: phonons (photons) are generated because the initial zero-point state |0 > is not zeroth 
state for the time-dependent destruction operators bj] there are phonons with frequencies Qj in |0 > at the time 
moment t ; the number of phonons equals ^ Nj(t), where Nj(t) =< Q\bf (t)b j(t)\Q >. 



(6) 



jj "jWi wucxc iijy^j . V\Vj 

Energy, which is generated in phonon system at the time moment t, equal: 

E ph (t) = 2 «Mt)(iV>(t) + 1 -)-J2^ = lj2 «r W| [ A(r) ' dr 



(fast oscillating terms are neglected). This energy comes from the local mode: E p h{t) — Ei(0) — Ei(t); Ei(t) ~ 
w ; 2 (3o(i)/2, Qo(t) is the mode amplitude, u>i is its frequency. This relation between the energies E p h and £7 gives an 
equation for Ei(t). Differentiating (||) and taking into account that dE P h = —dEi, one obtains Ei(t) » —~f(t)Ei(t), 
where 

. , Tthw 2 '^i^ 1 f UJk du> p(u>) piuji — ui) r ^, „ , .. , r f\ 

AuJl k^i J^-i |1 - w 2 Ei(t)G(u - cj;)/(2wf)| 2 L J 
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stands for the relaxation rate at the time moment t (ujq — 0, w m +i = wj/2), p(co) — — 3(G(w)) is the phonon de: 
of states (for uj > 0), G(lo — w\) = G(lu)G(lu — uj{) is the two-phonon Green's function, expression for Tk is given inE 
The second term in square brackets in (|^) is important if the initial energy is close to critical point. When deriving 
expression it was taken into account that for E\ w Ek 

l-Q{E k -EM)e-^ 
lk E^E® ■ (8) 



III. ODD LOCAL MODE IN A CUBIC LATTICE 



To illustrate the theory we calculate the relaxation of the local mode associated with a light substitutional impurity 
atom in a simple cubic lattice. Within the approximation of the nearest neighbours interaction the potential operator 
has the form 

* '=£££ ^(^.r. o) 

a n m— 1 

where a — x, y, z are the directions of crystal axes, n — (n x ,n y ,n z ) is the vector of the lattice sites, ri a is the vector 
of the site nearest to n in a direction, 

Rn a = y/(Ron a + r ai i J 2 + fj£ - ^ (10) 

is the operator of distance between the nearest neighbours in the a-direction, i?on Q is the distance between the 
nearest neigbour sites, fpft a — qpft — qpn a , qp is the /3-component of the displacement vector q^ of the atom n, 
= ^xn + ^yfi + ^zri ■ By expanding V in the power series of displacement operators one gets 

V = \Y. ( v ^A + Vinjlnc + VmJnJ«f( a + V^J 3 aHa ) + ..., (11) 

a,n 

where 

m 

^ q =e^ 3 4:v-2), 

= J £ ^; 3 4: } (^ - 2)(m - 4) . (12) 

m 

The potential considered does not take account of the covalent interaction which leads to the chemical bonding. 
This (covalent) interaction can, however, be easily included in calculations by introducing additional terms of the 

type v k a and v k a ; 

We make calculations for a local mode associated with a vibration of a light substitutional impurity atom situated 
at the origin of our reference frame . In this case the solutions of classical equations of motion, corresponding to the 
local mode, satisfy the conditions for the mode localized on the impurity atom at the site n — 0. As a consequence, 
all mathematical calculations are simplified, since one can neglect contribitions of other atoms to the classical motion. 
This allows us to neglect also the variations of the frequency of the local mode and of the constants of anharmonic 
interaction which are caused by the variation of the amplitude of the local mode. 

Altogether there are 21 operators of Cartesian coordinates (three for each of 7 atoms includung the central atom 
and 6 its nearest neighbours) which contribute to the relaxation of the local mode in this approximation. Therefore, 
the interaction Hamiltonian is 



21 

H' = Ql COS Ullt E V 3nn'q n q n ' 



(13) 



The matrix {V3 nn '} is determined by the coefficients of the quadratic operators of cubic anharmonicity in (|ll]) with 

f/3fi a =Cp + q f3n - q 0fiai , (14) 

where Cp are Cartesian amplitudes of the local mode; n denotes the numbers of the central atom and its nearest 
neighbours and the Cartesian components of their coordinates. We choose (see Figure 1): q\ — xq, q 2 — Vo, Q3 — Zo, 
94 = x lx , q 5 = y lx , q 6 = z la> , q 7 = X-i x , q 8 = y-l„, <?9 = Z-l m , <7io = U y , Qu = gi2 = h y , gi3 = V -i y , Qu = X-i y , 

Qi5 = z-i v , Qw = h z , qn = 9is = Vu, Qw = z~u, Q20 = X-i z , hi = V-u\ C = \JCx + Cy+ Cf- 

The matrix w — {V3„„'}/V3o for three directions in the crystal lattice is given below (V30 = Vso a ; k — 1 + V^ 0a /V30 ) . 

The operator H' with the given matrices w can be diagonalized analitycally. Since the corresponding final analitycal 
expressions for the eigenvalues are complicated, we restrict ourselves in giving their non-zero values for some chosen 
parameters k (see Table 1 and Table 2). 



/ 





-Mi 


Mi 


-Ma 


M 2 


-M 3 


M 3 \ 


-Mi 


Mi 



















Mi 





-Mi 
















-M 2 T 








M 4 













M 2 T 











-Mi 










-M 3 T 














M 4 





V 


Mj 

















-Mi) 



[100], 



,(2) - 



1 

71 



( 


-M s 


M 5 


-M 6 


M 6 


-M 7 


M 7 \ 


-Mi 


M s 

















Mi 





-M 5 














—Mq 








M 5 











M 6 T 











-M 5 








-Mj 














M 8 





\ Mj 

















-M 8 / 



[1101 



w 



(3) = _!_ 





-M 9 
M 9 
-M? 
Mi T 
-Mj 
V Mj 



-M 9 
M 9 













M 9 


-Mio 


Mio 


-Mi 














-M 9 














M 9 














-M 9 














M 9 















Mn \ 






-M 9 7 



[1111 














( l 







Mi = 





1 


:) 


, M 2 = 





1 


: ) 




Vo 














0/ 





/I 









^0 1 




M 3 = 








: ) 


, M 4 = 


1 


: ) 




Vo 


t 


0/ 




Vo 


0/ 



K 1 



1 1 



M 5 = 





/ 1 \ 






( 


1 

1 


1 \ 

: 


M 7 = 


1 




= 


i 









Vo i J 






Vi 





o) 




1 K 1 1 \ 






/ 1 


i 
i 


A \ 


M 9 = 


110 


, M w 


= 




i 






V 1 1/ 






Vi 





J 






I 1 1 












Mn = 


1 


1 









\k 1 1/ 

Table 1. Nonzeroth eigenvalues of the matrices w, k = 5. 



[100] 


[110] 


[111] 


±8.8175 


±7.4270 


±6.9176 


±2.2361 


±5.4146 


±4.6212 


±1.5003 


±2.1830 


±0.9218 




±1.2256 


±0.5774 




±0.8570 


±0.3825 




±0.8259 


±0.3586 




±0.7071 






±0.5813 





Table 2. Nonzeroth eigenvalues of the matrices u>, k = 1.5. 



[100] 


[110] 


[111] 


±3.1948 


±3.8230 


±3.9880 


±2.2361 


±2.1830 


±1.7728 


±1.2422 


±1.6619 


±1.0431 




±1.1114 


±0.5774 




±0.8570 


±0.1377 




±0.7071 


±0.1106 




±0.3364 






±0.1557 





To proceed further we must know the functions of the local densities of states p m (uj). To find these functions 
one should specify the parameters for our lattice and the substitutional impurity atom. For qualitative estimations, 
however, it is possible to involve suitable standard expressions for density functions which in many cases give good 
approximation. 

For qualitative estimations we use the following density functions: 

Pm(u) = \amPa(u) + \ om p (uj) , (15) 

where 

Pa (uj) = (3lu/ujI)^T^ (16) 
is the density function for acoustic phonons and 



Po{uj) = (8/TT(uj 2 - LUol) 2 )\/( UJ o2 - - LJ ol ) , UJ ol < UJ < UJ o2 . 



(17) 



is the density function for optical phonons, X am and A om stand for the corresponding contribution of the phonon 
bands to the m - th density of state, A am + X om — 1, u) a , uj \ and tu 2 are the limiting phonon frequencies for acoustic 
and optical bands respectively. 

) and ( |l7| ) make a correct account of the Van Hove singularities in the density of states of phonons 
in the crystal lattice. We should underline, that the two-phonon Green's functions G m (uj) which are relevant to our 
process, are less sensitive to details of local dynamics than the one-phonon Green's function G m {u>). 

In Figures 2-6 the critical behaviour of the rate of relaxation of a strongly excited local mode is shown. To simplify 
calculations we chose a system of units where h = k = 1. On the plots time is measured in units w" 1 , energy is given 
in units of vibrational quanta u> a , K is a parameter of interaction. 

The time evolution of the local vibration is illustrated in Figures 2, 4. The corresponding rates of the decay are 
shown in Figures 3, 5. As can be seen, the process of relaxation is determined by the set of parameters for interaction 
coefficients and by the initial energy of the excitation of vibration. Each relaxation jump leads to the release of 
about 50% of local mode energy in several periods of vibration. Depending on the initial conditions, the most stable 
vibrations can be either in [100] or in [111] direction. This effect of anisotropic relaxation is specific for the mechanism 
described. 

In Figures 2, 3 the initial energy (Ei(Q) — 700) and the critical values for directions [110] and [111] are close. 
Therefore, the vibration in [100] is more stable. The situation is just the opposite in Figures 4, 5 {Ei(0) — 1800), the 
lifetime of the vibrations in [100] and [110] is much shorter than that of the vibration in [111] direction. 

We carried out a set of calculations for different temperatures. In Figure 6 a comparative behaviour of relaxation 
for two temperatures is shown. Generally, for higher temperatures the damping of the local vibration is faster. 

Finally we should pay attention to the following features of the described process of generation of phonon bursts 
by strong local vibration: 

1. A remarkable time delay between the excitation of the local vibration moment and the phonon burst, generated 
by this vibration; this delay may be of hundreds or thousands of vibrational periods, see Figure 6 - relaxation 
in direction [111]. 

2. Quasimonochromatic spectrum of phonons, generated in the burst. 

3. Dependence on the direction of vibration. Separate calculations should be carried out for each type of the 
crystall lattice with a correct account of all possible contributions from vibrations of atoms. 

These effects may be used for the experimental observation of the process. 
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IV. QUANTUM INSTABILITY OF INTRINSIC LOCAL MODES IN ANHARMONIC MONATOMIC 

CHAIN 

The case of the ILM in a monatomic one-dimensional lattice is of special interest due to the following features: 

• The energy Ei and the frequency uji of the ILM are correlated. Explosion-like releases of energy will take place 
if Ei exceeds critical energy. Monatomic one-dimensional chain is the simplest model to investigate whether this 
condition is fulfilled. 

• The one-site Green's functional has only the imaginary part for the range of frequencies \u>\ < u> m , where u) m is 
the top band frequency. Therefore, the two-phonon Green's function is real; two-phonon anharmonic damping 
causes instability of the ILM of the frequency u>x < 2uj m . 

The potential energy operator which includes linear, cubic and quartic terms has the form 

^EEt(^-^. ( 18 ) 



r 

n r=2 



where U r is the operator of the longitudinal displacement of the n-th atom from its equilibrium position, K r are 
harmonic (r = 2) and anharmonic (cubic: r = 3, quartic: r = 4) springs. The operators U n satisfy the following 
equations of motion: 



d 2 U n _y^K r 
dt 2 ~ ^ M 

r=2 



(U n+1 - UnY' 1 -(U n - Un^Y- 1 . (19) 



Here M is the mass of atoms. We suppose that the ILM w ith frequency uj < 2ui m is excited at the time t = 
at the site n = and its nearest neighbors (u> m — 2y/ K2/M). Anharmonic interactions are supposed to be weak, 

satisfying the condition <C Aq, where ~ y/h/2u) m M is the amplitude of zero vibrations, Aq ~ -J K2 / K± is 
the amplitude of the local mode. 

To take account of the intrinsic local mode we introduce the operators U n in the form 

U n = A n (t) + i n + q n , (20) 

where the "classical" displacements A n are supposed to be nearly periodic functions A n {t) ~ a n {t) coswji, a n (t) and 
in(t) are slowly time-dependent amplitudes and shifts satisfying the equations 



oj; a n — K 2 (a n+1 — a n ) + 2K 3 (a n+ ii n+ i — a„£ 

•»n+l 



+ 3K 4 (^+i/4 + a„+i^ + i " a 3 J '4 - a„&) 



= K2(Cn+i - Cn) + K 3 (a? n+1 /2 + ? n+l - a 2 n /2 - + 
+ 3^4(^+1^+1/2 + S+i ~ a 2 n i/2 - , (21) 

where we neglected the contribution of higher harmonics ~ cos3w;t, cos5cj;t, . . and introduced for convenience 
new variables: a n = a n - a„_i, f„ = £„ - ^ n _i. 

Slow variations of a n and which come from quantum effects, are described by time dependent operators q n {t), 
which satisfy the equations 

= W n {t){q n+1 - q n ) - W n _i (*)(«„ - . (22) 

Here W„(t) = + v n 2 ^ + 2ui 3 ' ) costoit, where v n 2 ^ — 2K 3 i n + 3^(5^/2 + 1^) and wi 3 ^ = K 3 a n + ^K^anin determine 
the change of springs caused by the ILM. The terms ~ K 3 (q n — g„_i) 2 , K±A n {q n — q n -i) 2 and K^{q n — q n ~i) 3 are 
supposed to be small and they were neglected. Taking account of these terms is important if uji > 2uj mi when the 
two-phonon decay under consideration is forbidden by the energy conservation law. 

The Green's function of a perturbed lattice G m is expressed via the Green's function of a perfect lattice G™ (w): 

The rate of decay of the ILM is determined by ^ with p(u>) = 2tt~ 1 ^s(G(lo)). The properties of the ILM depend 
essentially on the parameter K 2 / K2K4. According to our estimations, for realistic systems it has the value betweem 
1.2 and 2.1. Relaxation jumps take place if the denominator in (]?]) approaches zero (if \v m G(uji/2)\ — 1). 

Further analitycal and numerical calculations based on the real parameters of interaction have shown that the 
condition \v m G(u>i/2)\ = 1 is fulfilled for the cases when both the central bonds and the next to central bonds are 
accounted. Therefore we may conclude that intrinsic local modes on perfect anharmonic chain are not stable. For 
the typical interaction parameters derived from rare gas solids potentials this conclusion is valid for the frequencies 
of the ILM lying on the interval 1.2w m . . . 1.7cj m . 



V. CONCLUSION 



We showed that quantum (and thermal) effects lead to the creation of extra channels in the decay of a strong 
local vibration (intrinsic or associated with an impurity) with the time-dependence being strongly nonexponetial. 
At very large energies the rate of the decay is relatively low. It increases during the course of the relaxation, at 
critical points it gets extremely high values. Then the rate decreases and after that increases again until the energy 
reaches next critical point. The energy of the local mode drops in a short time at each critical point and a burst 
of quasimonochromatic phonons is generated. The relaxation rate is higher for higher temperatures, although the 
values of critical energies do not depend on temperature. We also showed that the law of relaxation depends strongly 
on the directions of vibration and its behaviour is different for different initial energies. Strong dependence of the 
relaxation on the directions of vibration obtained for cubic crystals is remarkable, since in the case of the existence of 
chemical bonds the most stable are vibrations in the direction of the bonds. This effect may have an important value 
for chemical reactions and for the mechanism of defects formation in solids 



It follows from formulae (|7]) and (||) that the damping rate 7 of the local mode is strongly enhanced (diverges in 
our approximation as ~ \t — ife| _1 ^ 2 ) if the mode energy approaches (at t = tk) one of the critical energies E}.. This 
enhancement of 7 is associated with the generation of quasimonochromatic phonons. These phonons are emitted in 
pairs: one phonon with the frequency Uk and another with the frequency {u>i — u>k)- As a result the decay of the 
strongly excited local mode is highly non-exponential: it has step-wise jumps near critical energies. 

At a final stage of decay the relaxation becomes exponential, i.e. the rate of decay is constant. In weak coupling 
limit the theory presented also gives the constant rate of decay. This is similar to the result given by perturbation 
theorycj and to the numerical simulations!] of classical anharmonic chain. 

For experimental observation of the relaxation jumps both, quasimonochromatic spectrum of generated phonons 
and the time delay between creation of the strong vibration and the jumps can be used. Strong local vibrations can 
be excited e.g. by light or by high energy particles. Excitation by light takes place e.g. in luminescence centers in 
ionic crystals with small quantum yield of emission; here, due to the nonradiative transitions, the electronic excitation 
is abruptly transformed to the local vibrations. High energy particles create in crystals focusons which, in turn, 
efficiently pass their energy to impurities by excitation of their vibrations. 
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Figure captions 



Fig. 1 Displacements which contribute to the relaxation of the central atom. The central (impurity) atom and its six 
nearest neighbours are shown. 

Fig. 2 Relaxation of the energy of the local mode in directions [100] - thick line, [110] - thin line, [111] - dotted line. 
k = 5, Ei(0) = 700. 

Fig. 3 Rate of the relaxation of the local mode in directions [100] - thin line, [110] - thick line, [111] - dotted line. 
k = 5, 25/(0) = 700. 

Fig. 4 Relaxation of the energy of the local mode in directions [100] - thick line, [110] - thin line, [111] - dotted line. 
k = 5, 25/(0) = 1800. 

Fig. 5 Rate of the relaxation of the local mode in directions [100] - thick line, [110] - thin line, [111] - dotted line. 
k = 5, 25/(0) = 1800. 

Fig. 6 Relaxation of the energy of the local mode in direction [111], 
T = - thick line, T=0.85 - thin dotted line, k = 5, 25/(0) = 700. 
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